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Abstract 

We use the expectation of the range of an arithmetic Brownian motion and the 
method of moments on the daily high, low, opening and closing prices to estimate the 
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volatility is used to calculate Black-Scholes prices for European options, and a trading 
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1 Introduction 

This article is a modified version of what has been studied in the Ph.D. thesis of Kone 
(1996). It concerns the application of the method of moments to range-based volatility 
estimation using daily high, low, opening and closing stock prices. Aiming to estimate 
volatility and not to measure it, we assume a Black-Scholes framework with constant 
volatility and use daily data to achieve it (see Rogers and Zhou (2008) for further 
motivation for this choice). 

Incidental to this is the derivation of the density and expectation of the range of an 
arithmetic Brownian motion. Subsequent to this thesis, portions of it have been studied 
for different purposes (see, for instance, Sutrick et al (1997) for the use of the density 
of the range of an arithmetic Brownian motion in the do-nothing option, or Magdon- 
Ismail et al (2000, 2004) for the use of the expectation of the range of an arithmetic 
Brownian motion in different contexts). In particular, expressing the density of the 
range in the context of Sutrick et al (1997) corrects their expression. 

The literature on range-based volatility estimation includes classic work by Gar- 
man and Klass (1980), Parkinson (1980), Rogers and Satchell (1991) and Rogers et al 
(1994), whose estimators are reviewed in Yang and Zhang (2000). Of these, the latter 
paper is most related to the current one because it considers after-hours price jumps 
in addition to drift. However, the methods presented here are differenlH, and perhaps 
more practical (see the remarks of Chan and Lien (2003) on the empirical availability 
of certain parameters in Yang and Zhang (2000)). 

Starting from the joint density of the running maximum and the current value of 
an arithmetic Brownian motion, the density of their difference (referred to as half- 
range) is obtained, and its expectation computed. This allows the computation of the 
expectation of the full range (defined as maximum minus minimum), which will then 
be used in the method of moments for intra-day volatility estimation. 

After-hours arrival of information results in price jumps at the opening, and we 
model this as a virtual trading day which is unobservable, but which, when succeeding 
the trading day, gives on average the complete statistical representation of one day. 

Black-Scholes option prices are computed using the parameter estimates, and when 
they differ the most from the observed market prices a profit is made by an appropri- 
ately devised trading strategy. 

The paper is organized as follows. In Section 2 we derive the expectation and 
the density function of an arithmetic Brownian motion. In Section 3 the method of 
moments is used to estimate the parameters of the stock price based on daily high, 
low, opening and closing data. The estimated parameters are then used in Section 4 
to price European options on the stock, which are then compared to market prices to 
identify instances of flagrant differences. The effect of the mispricing is estimated by 
computing the profit to be made in these opportunities. We conclude in Section 5 with 
some comments on the efficiency of the method of moments in volatility estimation. 

^see Remarks ITT] and fT2\ 



2 



C. Buescu, M. Taksar and FJ. Kone Volatility estimation using high-low-open-close 



2 The range of an arithmetic Brownian motion: expec- 
tation and density 

Let {fi, T , P} be a probability space endowed with a filtration { J-f }t>o, and let { W^t}t>o 
be a one-dimensional standard Brownian motion adapted to {Tt\t>Q- For t > let Xt 
denote a standard arithmetic Brownian motion with drift /i and volatility a > 0: 

Xt = p.t^aWu Xo = 0, (2.1) 

and let {Aft}t>o, \m{\t>Q and {Rt\t>Q denote its running maximum, minimum, and 
range, respectively: 

Mt := sup X„ mt = inf X„ Rt := Mt - rrit. (2.2) 

Q<s<t Q<s<t 



First we derive the expectation E[Rt] of the range of the arithmetic Brownian mo- 
tion Xt. This is achieved by computing the density and expectation of the half-range 
Mt — Xt from the joint density of Xt and Mt. 

Lemma 2.1. The joint density function of an arithmetic Brownian motion and its run- 
ning maximum can be expressed as: 



P{Xt e da, Mt G db) = ^^^L=^ exp { -^^^-^ + 4« " 1^ da db. (2.3) 



Proof. The proof is standard. Using the martingale = exp{PV(/i/cr— tyU^/(2cr^)}, 

Girsanov's change of measure defines a new probability measure P for any measur- 
able set A by P{A) = E[Zt{Y)lA]. Theorem 3.2.2 of Karatzas and Shreve (1988) 
with Yt = jji/a and Nt = aWt gives that Nt = aWt — yut is a local martingale. The 
process Wt = Wt — t^/ a is a Brownian motion under the new probability measure P. 
Equivalently, aWt = jjt + aWt is a Brownian motion with drift under P, and we can 
write: 

P{Xt < a, Mt <b) = P{aWt < a, aMt < b) 

= J exp (^x - ^^t) PiaWt G dx, aMt < 6) .(2.4) 

An application of the reflection principle gives: 

P{aWt < a, aMt < b) = p(Wt<-^- p(Wt < -, aMt > b^ 



^1 \ 1 ^^'^b-a 
$ — ^ - 1 + $ 



ay 
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where we denote by 0(-) and $(■) the standard normal density and cumulative distri- 
bution functions, respectively. 

Differentiating the formula above with respect to a gives: 

PiaWt eda,aMt<b) = ^{M^)- ) 1 . (2.5) 




Replacing (12.51) in (12.41) and differentiating first with respect to a gives: 

P(X, e da. M, < 6) = ^ exp (i^a - i^*) (^(-^) - 

and differentiating then with respect to b gives: 

P{Xt G da, Mt G db) = ^exp (-^a - li^t) ( - -^)^'('^tL^) da db 
2 \ //i l/i^ \ 1 / l/26-a\2\ / In 26-a , „ 

^ j - 2 7!^ - 2 (ttt) K - 2 J 'ttt 

2(26 -a) r (26-a)2 /i 1 /iM , „ 
^ ^exp<^-^ + -^a- -^t} dadb. 



v^a3 ^ i 2tfT2 a2 2 a2 

The joint density /xt,Aft (a, 6) is given by the term multiplying da db above. □ 

Remark 2.1. This is one of those results that seemed to be always at hand (it can be 
obtained from equation (1.8.8) of Harrison (1985)), but never derived. Note the typo 
in Yang and Zhang (2000), whose expression (Bl) has a plus for the first fraction in the 
exponential. For a = 1 this was used in Example E5 ofKaratzas and Shreve (1998) in 
relationship to Clark's formula to obtain explicitly the hedging portfolio. 

The density of the half-range Mf — Xt can be obtained using a standard two- 
dimensional transformation of the above joint density. 

Lemma 2.2. The density of the half-range Mt — Xt is given by: 

Proof. For Yi = Xi + X2 and Y2 = Xi — X2 the joint density is /yi,y2(z/i; Z/2) = 
\fxuX2 . Taking Xi = Mt and X2 = Xt and using Lemma Ogives: 

r I N 1. (y\^V2 v\-y2\ 1 . (yi-y2 yi + yi 

jYuYiiyi, y2) - 7;JMt,Xt [ 7^ , 7^ ) - o JXt,Mt 



2 2/2 ""^^'^"^ V 2 2 

1 2yi + 2y2 - yi + y2 [ (2yi + 2y2 - vi + y2f yi-y2 1 /^^ 

2 V2^a^ ''''^l Sta^ 2 2a^ 

_ 1 yi + 3y2 f jyi + 8^2)^ l^ yi-y2 l/j^ 

2 x/2^a3 ''''^ I 8t(T2 +(t2 2 2 

1 yi + 32/2 / 1 (yi + 3^/2 - 2/it)2 u 

= 2 7i^'"n"2 — — 
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Note that Mt >0 implies Yi > —Y2, thus the marginal density of Y2 = Mt — Xt is: 

/oo 
fYuY2iyi,y2)dyi. 
■y2 

A change of variable z = {yi + 3?/2 — 2/it) / {2ay/i) gives: 
z > zq := — — dyi = 2a\/tdz, 

therefore: 



' 20 



This can be rewritten as: 



e-# exp {-2i^!,4 + 2i^*(-.„) exp {-2^!,4 



cr2 V cr^t 

2 r (/it - c) 

+ = exp ' ^ 



■|exp(-2^c)jdc, (2.7) 



(7V^ I 2ta2 
and the result follows. □ 
Proposition 2.1. The expectation of the half-range is given by: 



5 



C. Buescu, M. Taksar and FJ. Kone Volatility estimation using high-low-open-close 



Proof. A simple calculation yields: 




A change of variable c = za\^ — fxt gives the result. □ 

Consider now E[Xt — m^. For each path of the Brownian motion Xt with drift /i 
consider a symmetric path of a Brownian motion Xt having drift —yi,. Then Xt — vtit = 
— {Xt — Mt) and E[Xt — mt] can be calculated using the equation (12.81) with /i replaced 
by — /i. Whereas the formula for the expectation of the range follows. 

Theorem 2.1. The expectation of the range of the arithmetic Brownian motion Xt 
defined in ( 12. il) is given by: 




(2.9) 



Let us denote this expected range function by ER{ji, cr, t). On closer inspection 
this can be further simplified as a function of just two quantities: 

E[Rt] = ER{f^,a,t) = h(-^—), (2.10) 

where the function h is defined by: 

h{x, y) := + l)(2$(x) - 1) + exp ( - | y. (2.11) 

Note that ER{fi, a, t) = ER{fit, a\ft, 1) as it should (the range over a time interval 
(0, f) of an arithmetic Brownian motion with parameters /i and a is the same as that 
over (0, 1) when the parameters change to /it and a\Jt'). 

In the remainder of this section we derive the density of the range Rt of the arith- 
metic Brownian motion Xt. This is achieved by the use of the joint density of the 
minimum and the maximum of X^, a result with its own merit, that we could not find 
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published prior to Kone (1996) (Borodin and Salminen (1996, 1.15.4) published in 
the same year the joint cumulative distribution function only in terms of some definite 
integrals). 

A version of this result is used in Sutrick et al (1997) for the same purpose, but 
there it seems to have incorporated an error. 

To obtain the joint density F(a, h) of the maximum and the minimum we start with 
a lemma. 



Lemma 2.3. We can write: 



F{a,b) 



h{a, b, x) exp ( -^x — -^t ] dx. 



2 a' 



where 

h{a, b, x) 

hi{a, b, x) 



hi{a, b, x) — h2{a, b, x), 



(2.12) 



(2.13) 



h2{a,b,x) 



fc=— oo 

X exp 

oo 

E 

fc=— oo 



2k{2k - 2) 
aH^/2Trt 



[2k{b -a)-2b + xf 



[2k{b - a) - 2b + x] 



2ta^ 
[2k{b - a) - x]^' 



exp 



[2k{b -a)-x] 



Proof. Using the change of measure of the proof of Lemma 12. II and Girsanov's theo- 
rem we have: 



P{a <mt < Mt <b) = J pt{x; a, b) exp (^^x - 
which gives F(a, b) via: 

d'^P{a <mt < Mt <b) 



(2.14) 



F{a,b) 



da db 



Leibniz rule of differentiation: 



dz 



f{x, z) dx 



a{z) 



'^^^df{x,z) db , . .da 

dx + f{b[z), z)— - f{a{z), z) 



a{z) 



dz 



gives the partial derivative of (12.361) wrt b: 



, oo 



fc=— oo 



db 



2k{b — a) — X 



a 



db 



dz 



2k{b -a) -2b + x 
/i 1 /i^ 



dz 



(2.15) 



(2.16) 



X exp (^-^x - ^^t)dx. {Ill) 
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Differentiating wrt a this last equation gives: 

d^d) /2k{b-a)-x\ d^(p /2k{b - a) - 2b + x 



F{a,b) 



b ^ oo 



da db 



a^t ) da db\ o^t 



xexp rfx.(2.18) 



Direct calculation gives: 



(9V (2k{b-a) -X 



da db 



_d_ 

a^/t / da 
4P 



-2k)\2k{b-a) - x\ 



2nta' 



exp 



/27rt(T2 

[2k(b -a)- xf 



exp 



[2k{b -a)- xf 
[2k{b - a) - xf 



,(2.19) 



9V (2k{b~a) -2b + x 



da db 



d_ 

da 



-(2k - 2)\2k{b - a) - 26 + x] 



^k{k-\] 



exp 



exp 

\2k{b-a)-2b^xf 



Substituting (|2.19D - (|2.20| ) in (12.181 ) gives the resuh. 



[2k{b - a) - 2b + x] 
[2k{b-a) -2b + x 



(2.20) 



□ 



Proposition 2.2. The joint density function F{a, b) of Mt and rrit can be represented 
as: 



with 



F{a,b) = Fi{a,b) - F2{a,b) - Fs{a,b) + F^{a,b) 
- F^ia, b) + Fe{a, b) + Fr{a, b) - Fs{a, b), 



(2.21) 



FUa,b)= y ^^^^ ^ [(2k - l)b - 2ka + fit] 



k=—oo 



fx 



X exp <^ -^[2{k - 1)6 - 2ka] - 



[(2fc - 1)6 - 2ka - fxtf 



,(2.22) 



F2M= J-^ [2(A:-l)6-(2A:-l)a + /it] 



fc=— oo 



xexp<^ -^[2(A;-l)6-2A;a] 



[2(A;- 1)6- (2A;- l)a - fitf 
2t^2 



,(2.23) 
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A:=— oo 



xerf|P*-^"'-^'=°-'"|. (2.24) 

(TV 2t 



F^{a,b)= J2 



4k{k-l)fi'^ 



k=—oo 



2(7^ 



exp|-^[2(fc- l)&-2fca]} 



^^^^^2(k-i)b-(2k-i)a-,t^^ (2.25) 

(TV 2t 



4k 



F,{a,b)= J2 TTT^ [{2k + l)b-2ka + ^it] 



k=—oo 



X exp { —-^2k{b — a) — 



[(2fc + l)6-2A;a-/it]^ 
2t^2 



,(2.26) 



4A; 



F,{a,b)= V [2A;6-(2A;-l)a + /it] 

— tcy^y Ztxi 



k=—oo 



X exp <! -^2k(b - a) 



(2.21) 

[2kb- {2k-l)a- fit]"^' 
2to^2 



Fr{a,b) 



fc=— oo 



exp<{-^2A;(6-a) [> erf 



{2k + l)b - 2ka - fit 



aV2t 



,(2.28) 



^8(a,0)= 2^ ^:^-:^exp 



A:=— oo 



2cr4 



zi , ,T ^ I 2kb — (2k — l)a — ut , 

^ 2k{b - a) } erf | ^ j=J — ] ,(2.29) 



(T 



(Tv/2t 



where 



erf (x) = ^ / e * dt. 



(2.30) 



Proof. From Lemma |2J] we write F(a, 6) as the difference of two terms, which we 
denote Ii and I2: 



j hi{a, 6, x) exp {^^x - d^, 
■■ J /i2(a, b, x) exp (^-^x - 2 



(2.31) 
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In each Ji and I2 we combine the exponents and then use, respectively, a change of 
variable: 

X + 2k(b - a) -2b- ut x + 2k(b - a) - at 
z= ^ j= ^, z = ^ ^ ^ , (2.32) 

followed by an integration by parts for / ex'p{—z'^ /2)dz and replacement of 

by erf(x) via = 0.5eri{x/\/2) + 0.5. The resulting eight terms are then then 

denoted Fi{a, 6), z = 1, . . . , 8. □ 

We use the above expression to derive the density of the range. To make it suitable 
for comparison with that obtained by Sutrick et al (1997) in their Proposition 1, we 
change A; — )■ A; + 1 in the summations of Fi, F2, F^ and F4. 

Proposition 2.3. The density function f^^ (r) for the range of an arithmetic Brownian 
motion can be written as: 

-.00 _ 00 

/«,(r) = AkH{k) + ^ 5^ ^Kk + 1) (2.33) 

k=—oo k=—oo 

where 

I(k) = e-^(l + c2)(0(i^i-c)-20(J^o-c) + 0(ir„i-c)) 

+ e-^[{c^Ki -2c- c^)^{Ki - c) - 2{c^Ko -2c- c^)^{Ko - c) 
+ {c^K_i-2c- c^)^{K_i- c)], (2.34) 



and 



2/ifer 2/j(fc + l)r 

J(k) = e^{(l){K, + c)-(l){Ko + c))-e^^{(l){K2 + c)-(l){Ki + c)) 



~ 



+e-'^[-'-HK,-c)+'-HK 
-e-^^ ( - '-m2 - c) + - c] 

2nkr 

+e^($(/ii + c) - + c)) 

2n(fc + l)r 

-e^^mK2 + c)-^Ki + c)), (2.35) 



with 



{2k + 2)r {2k + l)r 2kr {2k - l)r fiVi 

K2 = ^ , Al = -= , Kq = -=, K„i = , C = . 

ay/t ayt ayt o\Jt cr 

Proof. After replacing k hy k+1 in Fi,i = 1, . . . , 4 of Proposition l2.2[ a two-dimensional 
transformation a = u — v,b = u, gives, via Jacobian, the density of the range and the 
running maximum. Its marginal density is the one we seek: 



fir) 



pr 

/ F{u — r,u)du. 
Jo 



Applying a change of variable and integration by parts gives the result. □ 
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Remark 2.2. This result corrects that of Sutrick et al (1997) where there appears to 
be a mistake in the computations. 

Remark 2.3. The probabilistic starting point for both Kone (1996) and Sutrick et al 
(1997) ispt{x]a,b) dx := P{a < nit < Mt < b,x < Xt < x + dx\Xo = 0). The 
former uses a result of Feller (1951 ) that can be traced to Levy ( 1948): 



^ oo 

Pt{x]a,b) = — -= V 



/2A;(6- a) 
V a^/i 



X 



2k{b-a) -2b + x 



a 



(2.36) 



while the latter uses a result of Billings ley ( 1968): 



^ oo 

Pt{x;a,b) = — -= 4)(^ 



'X + 2k{b - a) 



a 



2b-x + 2k{b-a) 



a 



(2.37) 



The probabilistic results (i2.36\) and A2.37\) are in fact equivalent, as one can be ob- 
tained from the other by appropriately replacing the summation index k with —k and 
(h(x) with (h(—x). 



3 The method of moments appHed to volatihty estima- 
tion using daily high, low, opening and closing prices 

In this section we apply Theorem 12.11 to the estimation of the drift and volatility pa- 
rameters of the stock price from market data on high, low, opening and closing prices. 

Definition 3.1. i) A trading day is the period elapsed between the opening and the 
closing bells of a calendar day. 

ii) A virtual trading day is the after-hours period beginning from the closing of one 
trading day and ending at the opening of the next trading day. 

Hi) A one-day period consists of one trading day followed by one virtual trading 
day. 

We assume that the stock price St has the usual geometric Brownian motion dy- 
namics: 

dS 

--^ = fis dt + a dWt, t>0. (3.1) 
St 

Then the log-stock price log St is the arithmetic Brownian motion Xt defined in (12.11) 
with drift coefficient ^ 

At = yu, - y . (3.2) 

Note that fig is the one-day period drift of the stock price St, while fi is the similar drift 
of the log-price Xt = log St, the volatility parameter a is the same for both St and Xt. 
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The market data used for parameter estimation is as follows: for each one-period 
day ?G{l,2,...,n}we denote by Si^i the opening price and by Hi and Lj the intra- 
day high and low prices, respectively (i.e. the high and low are observed only during 
the trading day, and not the virtual trading day - see Figure [I])- 

The after-hours arrival of information in the market determines a jump between the 
closing price of one trading day and the opening price of the next day. We model this 
jump by letting the same geometric Brownian motion St have an unobserved evolution 
during a virtual trading day. The length of this virtual trading day is assumed to be, on 
average, a fraction / of the unit length of the one-day period. 

Remark 3.1. This assumption follows Garman and Klass ( 1980) and Yang and Zhang 
(2000), except that they assume the after-hours trading day precedes the actual trading 
day. They call it the opening jump (from Ci_i to Oi), and assume it is modeled by a 
Poisson process. 

Thus for z G {1, 2, . . . , n} we have (see Figure[il): 

OPEN(^) = 0i = Si_i, CLOSE(«) = Q = 5^./, HIGH(z) = Hi, LOW(«) = L^, 

(3.3) 

where: 

Hi = sup St, Li= inf St. (3.4) 



The evolution of the price during the trading period i — l<t<i — fis given by: 

log St = logO^ + fx{t-i + l) + a{Wt - Wt-i+i), (3.5) 
and during the after-hours virtual trading period i — f < t < ihy: 

log St = log Ci + fiit-t + f) + aiWt - Wi.f). (3.6) 

Taking expectation in (13.51) when t {i — f) and in (13. 6|) when t i gives: 

i^flog^l =E\\og^] =/i(l-/), (3.7) 

E \og^ =i^riog-^l =^/. (3.8) 



Ci J L Si^f. 

Using Wt — Ws identically distributed to Wt^s, the trading day and virtual trading day 
variances are obtained, respectively, as: 



VAR 
VAR 



log 
log 



= VAR 


[log- 








Si^i\ 



a' a 



1 = VAR 


[log- 















(3.9) 
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Thus, we can write heuristically: 



VAR 



log 



VAR 



log 



O 



i+l 



VAR(trading day) + VAR(after hours) 



To estimate the variance over the trading day we use the method of moments. The 
range = log Hi — log Li of the arithmetic Brownian motion Xt = log St over 
the trading day [0,1 — /] was obtained in equation (12.101) : 



EiRi^f) = ERifx, a, 1 - /) = - /), a,/l~f, 1). 

In (|3.10l) we estimate E{Ri_f) using the daily range data: 



(3.10) 



i=l 



(3.11) 
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and/i(l -/)by (see (lO) ): 



..:4^1og^. (3.12) 

i=l 

This leads to the following equation to be solved for x, the estimate of a^l — f : 



h = h[^,^]. (3.13) 



2 

The squared of this solution gives an estimate Vi = x^ of the variance (volatility 
squared) corresponding to the trading day part of a one-day period. 

For the after-hours part of the one-day period we have two choices: Vq (centered 
approach) used in Yang and Zhang (2000), or (non-centered) used in Garman and 
Klass (1980). Using the former (i.e. the sample standard variance Vq), we obtain the 
estimate for the variance of the entire one-day period as: 

Vz■.= y^^y^1 (3.14) 

or, in annualized form, as: 

al ■= 2h2Vz. (3.15) 

Denoting by Vc the sample variance of \og{Ci/Oi) used in their estimator by Yang 
and Zhang (2000) : 

Vyz = + kVc + (1 - k)VRs. 

where Vrs is the estimator of Rogers and Satchell (1991) and Rogers, Satchell and 
Yoon (1994) and A; is a constant, we note the following. 

Remark 3.2. i) The term Vi replaces the linear combination of Vc and Vrs used by 
Yang and Zhang (2000) for the intra-day trading period, and it does not need estimat- 
ing the value ofk that achieves minimum variance. 

ii) Our estimator is a true range-based estimator (log-range to be precise since it 
uses \og{Hi/ Li)), unlike that of Yang and Zhang (2000). 

Hi) Our estimator Vz is independent of both the drift and the weight f of the after- 
hours information. 

Example 3.1. Consider the market data on the high, low, opening and closing prices 
for the IBM stock for the period from May 26, 2010 to June 18, 2010. For each of 
these days we consider the historical 3-month % estimates of ki and k2, and we solve 
the corresponding equation di.iil) . 

The solution is our estimate of the volatility — f corresponding to the trading 
day, and we present it in annualized form (i.e. multiplied by y252) in Figure^ We 
compare our estimate of the volatility corresponding to a one-day period with the one 
of Yang and Zhang (2000). On June 18, 2010 they are (Ja = 0.2781 (see ^3.15^ ) and 
0.2982 (annualized volatility corresponding to Vyz)- 



^For parameter estimation Hull (2006, p. 287) recommends using historical data of 90 to 180 days. 
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Figure 2: Estimated intra-day IBM volatility - May 26 to June 18, 2010 

4 European options: mispricing opportunities 

We use the resulting annualized volatility to compute the Black-Scholes prices of Eu- 
ropean options on the stock. We then seek those instances when the computed prices 
differ the most from the market prices, and devise trading strategies to take advantage 
of the price difference. 

We now use the volatility parameter estimated above to price European call options 
using the Black-Scholes formula: 

Ct = St - i^e-^C^-*) $(^2), (4.16) 

with 

^ \og{St/K) + {r + al/2){T-t) ^ 
aaVT - t 
d2 = di- aaVT - t. 

We devise a trading strategy to take advantage of the information differential be- 
tween our estimated prices and market prices. For simplicity we trade only in European 
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call options, and assume that at expiry there is a payment equal to the payoff so that 
no actual trading occurs in the underlying stock (naked trading). 

Having assumed a constant volatility there is no volatility smile and no stochastic 
volatilitjO so we restrict our analysis to European call options whose strike prices 
are relatively close to the stock price at the beginning of the period (preferably in 
the money), and whose expiry dates are up to three months (the parameters can be 
re-estimated later in view of new data). 

Example 4.1. Consider the market prices for the European call options on IBM for 
the period May 26, 2010 to June 18, 2010 with expiry dates June 18 and July 16, and 
strike prices K G {115, 120, 125, 130} (the stock price on May 26 was 125.91). We 
compare these market prices with the Black-Scholes prices calculated using A4.16\) . 
Here the inputs are the stock price, the volatility estimated in Example \3.1\ and r the 
value of the 1 -month US Treasury bill yield for the previous day (online\Treasury datJ^ 

Since the intra-day volatility of Example |3. II that we use in Black-Scholes formula 
does not include the effect of the after-hours evolution, we compensate by allowing 
our prices to differ by up to 10% from the bid-ask spread. Thus, we trade when our 
estimated call price falls outside the interval (0.9 x bid-price, 1.1 xask-price). 

There are two cases. If our price is lower, then we short-sell the option at the bid 
price and wait for the first day when the estimated price is no longer lower to buy back 
the option at the then ask price. If it expires and the call option is exercised then we 
buy the stock in the market and deliver it. 

If our price is higher, then we buy the option at ask price and wait for the first day 
when the price is no longer higher to sell it at the then bid price. If it reaches expiry 
date, then we exercise it. 

This trading strategy is summarized in Algorithm [T] for t between May 26, 2010 
and June 18, 2010 for European call options expiring at close June 18, 2010. The data 
is retrieved once a day, except on expiration date when it is retrieved several times a day 
(this can be implemented as an algorithmic trading strategy and deployed continuously 
without much effort, especially by those interested in technical trading). 

Remark 4.1. This strategy results in an overall loss of 0.68 (see TableU}. This is due 
mostly to one large loss induced by a large sudden move in the stock price on June 
10 (127.3 versus 123.9 the day before). That is because we use yesterday's intra-day 
volatility to trade in today 's world. 

Over a time horizon longer than a month the strategy can absorb such shocks in the 
stock prices, provided they are sparse. Alternatively, one can implement an additional 
stopping rule when the change in the stock price exceeds a pre-determined margin. 

A similar behaviour is exhibited when applying the same trading strategy to the Eu- 
ropean call option expiring at close July 16, 2010, but as the expiry date is longer than 
a couple of months the limitations of the assumptions of the model become apparent. 

Alternative approaches like stochastic volatility or econometric models (ARMA, GARCH etc) are 
not discussed here. 

^http://www.treasury.gov/resource-center/data-chart-center/interest-rates/Pages/default.aspx 



16 



C. Buescu, M. Taksar and FJ. Kone Volatility estimation using high-low-open-close 





K 


C(f) 




LA 0\X\> 






trade 


profit 


Mav 7ft 


1 10 


QO 


n 16 1 17") 


J til 


Mav 97 


CO 96 99") 


hi IV 

uu y 


17 

V/. 1 / 


Mav 98 


1 10 


57 


(0 75 78") 


J til 


Inn 9 

J Llll 


CO 67 


hi IV 

uu y 


05 


Jun 7 


125 


1.92 


(2.20,2.23) 


sell 


Jun 8 


(1.21,1.23) 


buy 


0.97 


Jun 7 


130 


0.29 


(0.42,0.44) 


sell 


Jun 8 


(0.15,0.17) 


buy 


0.25 


Jun 8 


120 


3.69 


(4.15,4.20) 


sell 


Jun 9 


(4.60,4.75) 


buy 


(0.6) 


Jun 9 


125 


1.09 


(1.30,1.38) 


sell 


Jun 10 


(3.05,3.15) 


buy 


(1.77) 


Jun 17 


130 


1.09 


(0.90,0.94) 


buy 


Jun 19 


(1.00,1.05) 


sell 


0.06 


Jun IS"" 


130 


0.60 


(0.48,0.51) 


buy 


Jun 18'' 


(0.63,0.69) 


sell 
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0.07 



at 12:27pm ''at 13:36pm at 15:58pm ''at 16:00pm if exercised 



Table 1: Trading in the call option expiring June 18, 2010 

(left: open a position, right: close position) 

5 Conclusions 

We have used the method of moments to estimate the volatility of the stock price and 
used this to identify arbitrage opportunities in the market of European options. As a 
by-product we have derived the density and expectation of the range of an arithmetic 
Brownian motion. 

In comparison to the estimate of Yang and Zhang (2000), our volatility estimate 
takes advantage of the actual range of the Brownian motion and perhaps does not 
overestimate as much. It is most useful for short expiration dates and for strike prices 
that are not far out. We believe it is an efficient alternative that can be easily computed 
and has a practical implementation. These traits recommend it to the attention of 
practioners in the field. 
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Algorithm 1: Trading strategy for a mispricing opportunity found at time t 



input : Parameters t, aa, r, St, bid(t) and ask(t) (European call prices) 
output: Profit of trading strategy 

1 profit^O; // initialize 

2 T ^ June 18, 2010; // expiry date 

3 compute CTc,; // Example 13 . II 

4 for 115 to 130 do 
C'(t) ^ BlackScholesCall (t,T,ii',cra,r, S'f) ; // (14 .161) 
if C < 0.9 X bid(t)then 

prof it ^ prof it -I- bid(t); // short-sell call 

while t < T and C'(t)<0.9xbid(t) do 

t ^ t-l-1 ; // wait 1 day 

compute (Tq; // Example 13 . II 

C{t) ^ BlackScholesCall (t,T, K,aa,r, St) / / (14 .161) 

end 

ift < Tthen 

prof it ^ prof it- ask(t) ; // buy back call 



else 



if call is exercised then 

profit ^ profit- (St -K); 
// buy stock and deliver for K 
end 



end 



se 



// buy call 



if C > l.lxask(t)then 

profit ^ profit- ask(t); 
while t < T and C'(t)>l.l xask(t) do 

t^t-l-1; // wait 1 day 

compute cTq; // Example 13 . II 

C{t) ^ BlackScholesCall {t,T, K, aa,r, St) ; 

end 

if t < Tthen 

I prof it ^ prof it -I- bid(t) ; // sell call 

else 

I prof it ^ prof it -I- (5*4 -K); // exercise call 
end 



end 
end 

35 end 

36 return prof it; 
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